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Наближене моделювання процесу кристалазаци 
при наявност! конвекци 


Рассматривается задача управления информационными процессами при автоматизации технологии 
тепловой обработки металла. В качестве источника информации исследуется математическая модель, 
основанная на пространственной задаче Стефана, с учетом конвективного движения в жидкой фазе. 
Существование гладкого решения проблемы доказывается посредством двух не слишком ограничительных 
допущений, касающихся исходных данных проблемы. Приближенное решение построено при помощи 
метода малого параметра. Построено приближенное решение задачи. 
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Розглядаеться задача управл1ння технолог!чним процесом теплово! обробки металу. В якост! джерела 
1нформацй досл1джуеться математична модель, базована на просторов!Й задач! Стефана, з урахуванням 
конвективного руху у рдиннЙ фаз!. Наближене рипення побудоване за допомогою методу малого 
параметра. Крум того, автором позначена проста умова, шо забезпечуе 1снування вдпов!дного стащонарного 
ришення 1 його стйкость, при цьому остання розумиеться в звичайному сенс!. Побудовано наближений 
розв’язок задачи. 

Ключов! слова: процес кристаллзацй, математична модель, завдання Стефана, радка фаза, 
конвекщя. 


Главной тенденцией в автоматизации является интенсивная разработка высоко- 
организованных систем управления технологическими процессами на базе совре- 
менных методов управления и средств вычислительной техники. 
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Стохастический характер функционирования производственных процессов, об- 
условленный случайным характером возмущающих воздействий и полезных сигналов, 
сложность математических моделей и критериев оптимизации приводит к существенным 
трудностям теоретического обоснования и решения практических проблем оптимальных 
автоматических и автоматизированных систем управления. На стадии проектирования 
АСУ ТП для уже существующих и строящихся печей возрастает важность обоснован- 
ности принимаемых инженерных решений. Разработка АСУ ТП должна опираться на 
исчерпывающие знания свойств процесса, изучаемого с точки зрения проблем автома- 
тизации управления. В этой связи рассматриваемая задача построения математической 
модели как источника информации о процессе является актуальной. 

Целью данной работы является обоснование математической модели как источ- 
ника информации для решения вопросов всестороннего анализа и управления инфор- 
мационными потоками при автоматизации технологических процессов нагрева металла в 
проходных печах прокатного производства непрерывной разливки стали на основе 
математического и имитационного моделирования, анализа статистических данных и 
теплофизических экспериментальных измерений. 

Работа посвящена изучению процессов кристаллизации двухкомпонентных сред 
в случае, когда распространение тепла связано не только с теплопроводностью, но и с 
конвективным переносом, присутствующим в жидкой фазе вещества. Рассматриваемая 
задача включает в себя как двухфазную задачу Стефана, так и начально-краевую задачу 
для системы Навье-Стокса, описывающую движение вязкой несжимаемой жидкости в 
нецилиндрической области. При изучении задачи учитывается скачек плотности 
вещества на границе раздела фаз. 

Для описания поля скоростей в зоне поступления перегретого металла исполь- 
зуется математическая модель затопленной струи вязкой жидкости, основанная на 
известном в теоретической гидродинамике точном решении нелинейной системы 
дифференциальных уравнений Навье-Стокса, а поскольку в затопленной струе пере- 
мешивание можно охарактеризовать так называемой свободной турбулентностью, 
характеризующейся одним числовым параметром — коэффициентом «кажущейся» или 
турбулентной вязкости Г,, полученные формулы можно интерпретировать как в 


ламинарном, так и в турбулентном приближениях. 


Постановка задачи 


Теплофизические процессы в кристаллизаторе, сопровождающиеся фазовыми 
переходами вещества, описываются математической моделью, в которой температура 
каждой из фаз удовлетворяет уравнению переноса тепла со своими теплофизическими 
коэффициентами. На границе раздела фаз обе температуры постоянны и равны тем- 
пературе фазового перехода (для химической однородной среды), а на заданных частях 
границы — стенках кристаллизатора, поддоне — поддерживается определенный режим 
(теплоотвод, теплоизоляция др.). Поверхность раздела фаз (фронт кристаллизации) 
является неизвестной, или «свободной», границей, и для ее определения дополнительно 
задается «условие Стефана», означающее, что тепловой поток через фронт кристаллизации 
в сторону твердой фазы равен тепловому потоку со стороны жидкой фазы плюс 
скрытая теплота фазового перехода. Жидкая фаза рассматриваемого процесса заслуживает 
специального исследования ввиду априорной возможности существования поля скорос- 
тей, вызывающего интенсивную теплопередачу путем конвекции. Усиленная циркуляция 
в расплавленной шлаковой ванночке была обнаружена в исследовании академика 
Б.Е. Патона и его сотрудников [1]. 

Цель состоит в изучении гидродинамических явлений в жидкой фазе. 
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1 Рассмотрим область © = 4х, „х ое ех фа В?) и через Г` и Г* 

обозначим следующие сферы: 
= аа. 2 2 уд + _ ле. 2 2 —_р2 
Г ={@, хх +) +; =Р УГ” = (<, х,,х,):х +, +: ЕЁ}. 

Далее, пусть Гу гладкая, связная поверхность без самопересечений, держащая 
внутри ©, которая разбивает ее на две подобласти О*'и О’, т.е. @=О* (]О` причем 
сфера Г’ лежит внутри ограниченной области, границей которой является Гу. 
Рассмотрим краевую задачу со свободной границей Г,. Требуется определить тройку 
(и (х),Г,) по следующим условиям: 

2. + +. + + . 
У?и* (х) =0,хеО*;и* (х) |= В* (5); 


и (х) ЕТ Ум” (х)- Ми] =0,хЕГ.. (2 


При этом, В*(х) и и* (х) — гладкие функции, а Г, принадлежит классу С” [2]. 
Затем введем в рассмотрение функцию и(х), заданную следующим образом 


и=и (>), при хЕО и и=и' (х) при хеО*. Тогда функцию и(х) можно найти из 

условия минимума функционала Г[(и,Г,)= ГУЧ ча на соответствующем 
[Ф) 

множестве А допустимых функций [2]. Это следует из формулы первой вариации 


интегрального функционала с неизвестной областью интегрирования [3]. 
Далее, удобно представить функционал / в сферических координатах: 


2лл КЮ 
_ 2 1 2 1 2 а 
Риго) = а ар В 5Ш 9афабар. (2) 


Пусть тройка (и*,Г,) является классическим решением задачи (1). Тогда эта 
тройка будет стационарной для функционала (2) множества К. 

Обратно, каждая стационарная тройка (и*(х),Гу) функционала на множестве №, 
где Гу — достаточно гладкая, связная поверхность, является решением задачи (1). 

Сформулированная задача (1) получается из задачи, изученной в [2], в случае 


Г =0, т.е. в случае бесконечно большой вязкости, Ве =0. 
Поэтому в дальнейшем под решением задачи (1) при Ке=0 будем понимать 


функции Го) = 0,м* (х) ии (х), заданные в ©*. 

Из условий (1) следует, что Г, — не что иное, как линия уровня функции 2х), то есть: 

Г ={хеЕО:и(х)=Ц. 
Далее, если предположить выполнение следующего условия: 
+(В*(х)-0> =, >0,хЕГ*, 

где =, — некоторая постоянная, тогда поверхность Г, лежит внутри области <> и 
представляет собой поверхность класса С”, не имеющую самопересечений и рас- 
полагающуюся относительно Г” и Г’ аналогично поверхности Г, (свободная по- 
верхность), изученной в [2]. Следовательно, рассматривая функцию и(х) в одной из 


областей Ф › И принимая во внимание лемму о нормальной производной, находим что 


==> 0.хеГ,, 
и 
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где и — нормаль к Г, направленная в сторону ОФ, , а = - некоторая постоянная. 
Отсюда, применяя теорему о неявной функции, следует, что Г, принадлежит классу 
С”, так как этому классу в некоторой окрестности Г, принадлежит гармоническая 
функция и(х). 

2 Минимум иная (2) на множестве А будем искать при помощи сумм: 

ве (в В ера Сир ьф.в) 

где у,($,0) -— ре функции. Неизвестные коэффициенты С, опре- 
деляют при помощи метода Ритца. Тогда поверхность Гу : р = р, (ф,0) определяется 
из уравнения и, (ф,0, р, (ф,0)) =1. 

При этом, необходимо учесть, что №и(> а: О, где =, — некоторая 
постоянная [3]. 

При малых / справедливо представление: 


и (Ф,0,10) 
Г, : р(ф,0,1) = ру(ф,0) - Ве 
Ю(ф ро (Ф У 


Здесь Ве — число Рейнольдса, а и! (ф,0,Г) — первое приближение исходной 
задачи, изученной в [2]. 
В частности для нулевого приближения и, (ф,0) из уравнения: 


+ 0(Ке), (ф,0) Е Г.. 


2 
и =” В _ В*) +(р’-"”)(В* - р^)С, =1 легко найти поверхность 


(Си) 8 < (В +" 5-5 Вбр а 2 (В” _ В*)+В*) -4С.В?т? 
р(Ф, 9) = == . 
При этом коэффициент С’ находим методом Ритца. Имеем: 
р? — 92 в й 
и = РВ —В,)=Р, 
В -г 


вы - + _ + + + 
ит (В, В) = ВЧ” +208? -)р=Е, (В В+) =Е,,В* = В*(ф,6,1). 


Таким образом, получим 


2лл К 
1 т 
Ци) = С В, +Е т ЕЙ, эт бЯфа ду. 
(и)= П. В РЕ т «404 
Построим уравнение Ритца: 
6) 2лл В 


=2 [ [С.Б + ЕДЕ? р’ эт бифавар = 0. 


0900г 


дс, 


Отсюда получим 
2х 


П 
00 
2лл В 2лл В 


с, =- [Е.Б р’ этбйавар | | [[Е? р’ этбйфавар 


00г 00 


Е.Е. р’ зт@афа бар, 


—_> 
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ИЛИ 
_ ИЕ 4’ +2(8В? — =? у: — В*) р’ за@Яфа@Яр 
[@ — 00хг 
| Е — 12 р]? р? за@ЯфАвЯр 


0900г 
На рис. 1 представлена поверхность Г,, которая расположена между сферами 


радиусов А и г при следующих значениях параметров: [= 200, = 0,8, 
-л/2<0<л/3,-л/2<ф<л/2,В* =3с05’0+с05° 6], В =-0,35[соз? 09 +соз? Ф]- 0.1, 


1 
и, (ф,0,1) ==» ф,6) +1. 


Рисунок 1 — Алгоритм построения поверхности Г, 


Предложенный алгоритм построения поверхности Г, позволяет исследовать 


численно эту поверхность в зависимости от основных параметров задачи. 


Литература 


1. Патон Б.Е. Избранные труды / Патон Б.Е. — Киев : Институт электросварки им. Е.О. Патона НАН 
Украины, 2008. — 893 с. 

2. Миненко А.С. Приближенный анализ многомерной конвективной задачи Стефана / А.С. Миненко, 
А.И. Шевченко // Допов1д! НАН Укратни. — 2010. — № 4. - С. 30-34. (фахове наукове видання 
Укра1ни, включене до затвердженого ВАК перелку). 

3. Миненко А.С. Методы исследования нелинейных математических моделей / А.С. Миненко, 
А.И. Шевченко. — Ки\в : Наукова думка, 2012. — 130 с. 


Г цегабага 


1. Раюп В.Е. З@есе4 ууогК$ / Раюп В.Е. — Клеу : Месёлс \Ме!Чте ш5Ийще патед Е.О. Раюп ое МА$ 
оРОКгаше, 2008. - 893 $. 

2. МшепКо А.5. Арргохипае апа|уз1$ оЁ ши-айпепз1опа| сопуесНоп оЁ З{еРап / А.5. МшепКо, 
А.Т. ЗВеусКепКо // Верог5 оЕ МАЗ оЁОКгаше. — 2010. — №4-5. 30-34. (ргое$$1юпа! }оигпа!$ ОКгате 
шс[а4деа ш Фе |184 оРарргоуеа НАС). 

3. МшепКо А.5. Мефо4$ Юг фе з4у оЁ попПпеаг тафетайса| то4е! / А.5. МштепКо, А.Т. ЗВеу- 
спепКо. — Клеу : МачКоуа ОиткКа, 2012. — 130 $. 


84 «Искусственный интеллект» 2013 №2 


Приближенное моделирование процесса кристаллизации при наличии конвекции 


ВЕЗОМЕ 


Т.А. Зурко 
Арргохтае Моаейиз ор ше СгуядаШ2айоп Ргосе55 


т Фе Ргезепсе о} СопуесПоп 
ТЬ$ рарег ехеп4 0 Ите-дерепепе сазе зоте гезаЙй оМашей Бу Фе аиог г 
збеаЧу-5еае еп ргоМет УИ сопуесНоп. Кеепте 1 Ще 1се-муаег зузет, уг 15 


аззитед ю Бе шсотргезЫе ап4 © обеу фе ЗюКез едиайоп =. — УЛу+У,р= Г (и) (у - 


(+) 


уе]осйу, у — КшетаНс у1$созИу, р -— ргезяше, Г — Баоуапсу Югсе, и” — умег 


{етрегаиге), \ВИе {етпрегацие и“ зайзНез 1е Веай соп4исНоп-сопуесНоп едцаноп ув 


{етрегаиге- дерепдепе Пеглорруз!са! ргорегйез. ТВе {етрегацие Не! и“ ш \е зоПа 
р8азе 15 соуеге Ъу Аазюп ошу. А Фе 1се-маег пиегсе, и“) =иГ) = 0 ап Ше 
З(еГап соп оп Во14$. ТВе зВете 1$ сотр! ее4 Бу ша апа Боипдагу соп41оп$. 

ТБе аиог ргезеп ап ех15${епсе Феогет Гог Фе сазе оЁйгее зрасе Читепз1опз. Те таш 
ЧЕЯсиу со15$1565 ш Фе ЁсЕ Фа ю шегрге! Фе ЗюКез едиаНоп ш Фе у’еаК зепзе, зоте 
шЮптайоп 1$ пееде4 оп Ше теслоп УТеге Ше 1етрегайе 1$ розшуе, умей 15 ш ит 
шИЯцепсе Бу Фе уеюсйу Не изе!. ТВе ргес1зе Ююгииайоп оЁ Фе ргоМет гедашез а 
{есиса|! сбое оп ЁпсНоп зрасез. Ех1\епсе оЁ зоайоп 1$ ргоуе Бу шнодисше а 
{етрегайте 4ереп4еп{ репаМу {епт ш Фе Яша ПНо\у едианоп ш ог4ег №0 4ейпе Бо Ше 


арргохипайпе {етрегацие и, ап арргохипайоп уеюсйу у, ш №е уБое доташ. Сот- 


распез$ агоитеп 5 аге изе4 {0 сеЁ а сопуегоепЕ зибзедиепсе, \Позе Пай 1$ зЗВо\т тю зОуе Фе 
опотпа] рго Мет. ТВе ацезвоп оРишачепез$ 1$ [ей ореп. 


Статья поступила в редакцию 09.04.2013. 
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